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. A direct method is developed to reconstruct the equation of state for high-density nuclear matter 

' from the relationship between any two properties of neutron stars, such as masses, radii, moments 

^jQ, of inertia, baryonic masses, binding energies, gravitational redshifts, and their combinations. 
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, The equation of state (EOS) for supranuclear-density matter with zero temperature is relevant to neutron stars. For 
^vq supranuclear-density matter, a number of possible scenarios, such as Kaon condensation, pion condensation, hyperonic 
matter, and strange quark matter, have been proposed. In the near future, lattice QCD is expected to predict the 
T-H • correct EOS for supranuclear-density matter. However, an accelerator-experimental test of the theoretical predictions 
for the EOS for cold high-density nuclear matter is extremely difficult. 
. In this light, neutron stars are very unique laboratories. Neutron stars are observed as pulsars. X-ray compact 
sources, soft 7-ray repeaters, and gravitational wave sources. Binary pulsar observations can determine masses of 
neutron stars with very high accuracy [Q. From observations of type-I X-ray bursts, kHz quasiperiodic oscillations 
(QPO's) of low mass X-ray binaries, and so on, some constraints on the radius of neutron stars have been obtained. 
T— I Actually, the value of the radius itself has been already measured (|] . Very recently, it has been proposed that the 
radius of neutron stars can be directly measured by the observation of gravitational waves from coalescing neutron star 
(NS) /black hole (BH) binaries [HJJ] . Some combinations of masses and spins of coalescing NS /NS and NS /BH binaries 
^ can be determined from gravitational wave observations in the inspiralling phase . Recently, it has been pointed out 
that the broad peaks at frequencies ^ 20 — 40 Hz of kHz QPO's can be interpreted as the precession frequency due to 
the Lense-Thirring effect, the relativistic frame-dragging effect due to the spin angular momentum This suggests 
the possibility of direct detection of the ratio of the moment of inertia to the mass of neutron stars. Observations 
' of absorption lines in X-ray spectra may reveal the gravitational redshift In some situations, from knowledge 

?• of progenitors of supernovae, we may infer the baryonic mass of the remnant neutron stars. In next decade, we will 
. ] have a good amount of data on the properties of neutron stars. 

. The structure of neutron stars has been investigated for a large set of EOS's. Arnett and Bowers Q examined 
?H ' systematically the mass-radius relation of neutron stars for various EOS's. Recently, Lattimer and Prakash |ic| ] 
examined systematically the relation of masses, radii, moments of inertia, and binding energies for a large number of 
modern EOS's. An alternative approach is a direct reconstruction of the EOS from the observed properties of neutron 
stars. Lindblom developed a method to reconstruct the EOS for nuclear matter from neutron-star masses and 
radii. In this report, we give a direct method to reconstruct the EOS from the relationship of any two properties of 
neutron stars. We adopt geometrized units such that c = G = 1. 

Before the method of reconstruction is derived, it is helpful to briefly review how observable quantities of neutron 
stars are obtained for a given EOS. For simplicity, we assume that neutron stars are static and spherically symmetric 
with no magnetic field. The line element is given by 

ds^ = -e'^'^'-'^dt^ + e^^^'Ur"^ + r\de'^ + sin^ Odcf^). (1) 

The structure of neutron stars is determined by the Tolman-Oppenheimer-Volkoff (TOV) equations given by 
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r(r — 2m) 

m' = Aner'^ (3) 

where P and e are the pressure and energy density of the matter, m is defined as 

l--=e- (5) 
r 

and the prime denotes the derivative with respect to r. For a given EOS, it is easy to integrate the TOV equations 
with the initial values P{0) = Pc and m(0) = 0. The radius R is the value of r such that P = at r — R. The mass 
M is obtained by AI — m{R). The moment of inertia / is calculated as 



6 t^{R) 



where uj{r) is obtained by integrating the following ordinary differential equation 
with the initial value lu{0) =/= and lu'{0) = 0. The baryonic mass msN and binding energy B are obtained by 



msN ^ / dre-^Anpr'^, (8) 
B = ms^ - M, (9) 



where p is the rest-mass density that is written as p = tobti using the baryon mass ruB and the baryon number 
density n. The gravitational redshift z at the surface is given by 

■ 2A/^ 



-1 (10) 

Then, we describe how to reconstruct the EOS from a complete observed data set. The problem is to determine 
the EOS p = p{P) for the domain Pq < P < Pu, where Pu is some cutoff pressure. We assume that we have a 
complete relationship f{Q^^\Q^'^^) = between observables Q'-^-' and Q'-^-' and that we have the correct EOS for 
P < Pq. We also assume that p is a nondecreasing function of P f\. We prepare N grid points in the domain [Pq, Pn] 
a.s Pq < Pi <■■■< Pi < Pi+i < ■ ■ ■ < Pn = Pu- Here, as a simple example, we take the grid points as 

log P, = log -Po + ^ (log Pn - log Po). (11) 

Suppose we have the "pointwise exact" EOS for P < Pi. This implies that if the TOV equations are integrated 

with the initial values P(0) = Pk and m(0) = then the obtained q]^^ and Q^^^ satisfy f {Q')^\ Q^^'' ) — *-* l^'^ 
< fc < z. Using some unknown number pi+i = p(Pi+i), we interpolate the EOS in the interval [Pi,Pi-\-i\. If it is 
taken into account that the EOS for nuclear matter is locally close to polytropic, then, the best interpolation is the 
logarithmically linear interpolation 

logP = logP, + r,(logp - logpO, (12) 
where Ti is the "local polytropic index" given by 

_ logP,+i - logP, 



logPi+i - logpi 



(13) 



Otherwise, the fluid is locally unstable. 
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Using this interpolation, e is obtained as functions of P S [Pi, Pi+i] as 



e = p{l + e), 



(14) 



where 



P= 17- 



K, 



1-r: 



(15) 
(16) 
(17) 



e is the specific internal energy and Eq. (|T^) is obtained from the first law of thermodynamics. Now that we have 
the EOS for P < P^+i parametrized by one parameter pi+i, we construct the neutron star by integrating the TOV 



equations with the initial values P{0) = Pi+i and m(0) = 0. Since the obtained qI^i and Q\^i are functions of 

Pi+i, we denote them as Q['^i{pi+i) and Q|^\(/9i+i). It is easy to numerically find the root p^+i (> pi) which gives 

f{Qi+iiPi+i),Qi+i{Pi+i)) — 0, by the bisection method or Ncwton-Raphson method with the aid of numerically 

estimated derivative df{Q^j^-^^(pi-^-l),Q\'^^{pi^l))/dpi+l. It is noted that {pi+i — pi) may not always be small even for 
small {Pi+i — Pi) because of the possible first order phase transitions. Thus, we find e^+i as 



l(2) 



K,, 



Si+l = 



) + e. 



(18) 



Then we obtain the pointwise exact EOS for P < Pi+i- 

Repeating these processes, we reconstruct the EOS for P < P^ — Pij. The reconstructed EOS is pointwise exact 
in the sense that the relationship /(Q*-^\ Q*-^-*) = is exactly satisfied by all neutron stars with the central pressure 
P € {Poj^ij''' ^W}- By taking the limit N oo, the reconstructed EOS becomes exact and in this limit the 
reconstruction does not depend on the interpolation method. 

In practice, it is very important how fast this scheme converges as N is increased. It is clear that each zth 
step reproduces the correct answer if the EOS in the interval [Pi,Pi^i] is a polytropic one. This implies that our 
extrapolation is linear in each interval with respect to logP and log p. For an estimate, we assume that logp is Taylor 
expandable with respect to log P as 



logp = log Pi 



^^1°^^^ (logP-logPO ' ' 



d? log p 



^dlogPy .'""^^ ' 2 V(dlogP)2 

The error produced in the interval [P^, Pi+i] is estimated for sufficiently large N as 

log p 



(logP-logP)' 



(error) 



,(dlogP)2 

Then the totally accumulated error is bounded from above as 



(logp+i -logp; 



N-l 

E 

i=0 



{error)j 



^ T7 (log Pjv - log Po ) sup 

^ IPo,Pn 



log p 



(dlogP)2 



(19) 



(20) 



(21) 



Therefore the convergence is as fast as 1/N. 

Next, we describe the reconstruction from a finite data set. In this case, we have only a set of N reliable data pairs 



{{q['\q?\ 



(22) 



The order of data pairs should be determined so that the smaller value of index i corresponds to a neutron star with 
the lower central pressure. The order of data pairs does not seem to be determined a priori except for the fact that 
M and msN are increasing functions of the central pressure for stable neutron stars. Unlike the previous case, we 
have to determine two unknown numbers Pi+i and pi+i such that 
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(P,+i,P,+i) = Q: 

Qi+liPi+lT Pi+l) = Q 




,(1) 



(23) 
(24) 



where p^+i) and Ql^i{Pi+i, pi+i) are the values of Q'^^^ and Q^^\ respectively, which are obtained by 

integrating the TOV equations with P(0) = P.+i and m(0) = with the EOS interpolated by Eq. (|3|)-(|l^ using the 
unknown numbers Pi+i and pi+i- This can be done numerically by the two-dimensional Newton- Raphson method 
with the aid of the numerically evaluated Jacobian 



Then we obtain the pointwise exact EOS for P < Pi+i. If a set of logP^ has roughly homogeneous distribution, the 
convergence is as fast as 1 /N as we have already seen. 

Here, we briefly discuss desired observational accuracy in order to draw some useful conclusions in discriminating 
between the possible EOS's, based on the figures given by Lattimer and Prakash The data pairs on masses and 
radii of neutron stars will be very powerful if the measurement is as accurate as AM ^ O.IMq and AR ^ 1 km. 
For the pairs I/MH^ and M/R, A{I/MR^) < 0.01 and A{M/R) < 0.01 would be desired. For the pairs B/M and 
M/R, A{B/M) < 0.01 and A{M / R) < 0.01 would be desired. However, the above estimate is quite naive, and more 
systematic and elaborate study is needed. The present accuracy for mass measurement from a radio binary pulsar 
is AAf ~ 0.0007Mq [Q, while that from an X-ray binary is AM ~ 0.4Mq For the measurement of radius, 

AR ~ 5.2 km is obtained if the interstellar absorption of X-ray emission and the distance to the object are fixed [|. 
More detailed X-ray spectral analyses and/or the space interferometric mission will dramatically improve the accuracy 
in determining the radii of neutron stars. Although moments of inertia of neutron stars have not been determined 
directly, the determination of them will be a complementary and potentially powerful technique. 

Finally, we compare the present method of reconstruction from a complete data set with that of Lindblom 's 
one. (1) The former applies to any two properties of neutron stars while the application of the latter is restricted 
to the mass-radius relation. (2) The former needs iteration to determine Pi+i for each zth step, while Pi+i and e^+i 
are determined explicitly and no iteration is needed in the latter. (3) In the former, the relation /(Q^^-*, Q*-^-*) — 
is satisfied at any grid point P E {Pq, Pi, • • • , P/v} for finite N, while it is not in the latter. (4) In the former, the 
obtained EOS is of the form P = P(p), which is conventional in nuclear physics, while in the latter it is of the form 



In summary, a complete data set of two properties of neutron stars is sufficient to reconstruct the EOS for cold 
nuclear matter. We have proposed a direct method to reconstruct the EOS. It is crucial to measure more than one 
property of a neutron star in order to put stringent constraints on the EOS. 
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P = P{e). 
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